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aids  the  elastic  deformation  of ^ much  more  strongly  on  the  basal 
slip  systems  than  on  either  prism  or  pyramidal  slip  systems.  The 
significance  of  these  interactions  stresses  on  the  initiation  of 
plastic  flow  is  considered.  The  effect  of  the  interface  phase  on 
the  elastic  interaction  stresses  has  been  considered  and  a  possib 
role  of  interface  stresses  on  elevated  temperature  creep  of  ^70^— 
alloys  has  been  suggested .  /  ^  \  >j 

By  use  of  a  NASTRAN  (18)  Computer  Program,  the  Finite  Element 
Method  (FEM)  has  been  employed  to  calculate  the  effect  of  particl 
size,  matrix  and  volume  fraction  on  the  stress-strain  relations 
of titanium  alloys.  It  was  found  that,  for  a  given  volume 
fraction,  the  calculated  stress-strain  curve  was  higher  for  a 
finer  particle  size  than  for  a  coarse  particle  size  within  the 
range  of  the  strains  considered.  This  behavior  was  seen  for  all 
the  different  volume  fraction  alloys  considered.  The  calculated 
stress-strain  curves  for  four  different  vol.  %  "(alloys  were 
compared  with  their  corresponding  experimental  curves  and  in 
general  good  agreement  was  found.  Whenever  there  were  discrepan 
they  were  discussed  by  comparing  the  morphology  of  the  mesh  used 
in  the  calculations  with  the  morphology  of  the  actual  materials. 
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PREFACE 


This  report  is  written  with  two  independent  parts,  which  are  en¬ 
tirely  self  contained. 

Part  I.  is  the  Role  of  Elastic  Interaction  Stresses  on  the  Onset  of  Plastic 
Flow  for  Oriented  two  Phase  Structures. 

Part  II.  is  the  Finite  Element  Method  (FEM)  Calculations  of  Stress 
Strain  Behavior  of  Alpha- Beta  Ti-Mn  Alloys. 
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PART -I:  THE  ROLE  OF  ELASTIC  INTERACTION  STRESSES  ON  THE  ONSET 
OF  PLASTIC  FLOW  FOR  ORIENTED  TWO  DUCTILE  PHASE  STRUCTURES 


ABSTRACT 

Elastic  interactions  at  the  Widmanstatten  or-3  titanium  alloy  interfaces 
that  arise  due  to  compatability  requirements  have  been  calculated  for  vari¬ 
ous  orientations  of  the  interface  with  respect  to  the  stress  axis.  It  is  shown 
that  maximum  interactions  are  found  when  the  stress  axis  lies  in  or  close 
to  the  interface  and  40°  away  from  [  0001 .  The  interaction  stresses 
could  be  of  the  order  of  30-35  percent  of  the  resolved  shear  stress  on  the 
basal  slip  systems.  2  aids  the  elastic  deformation  of  a  much  more  strongly 
on  the  basal  slip  systems  than  on  either  prism  or  pyramidal  slip  systems. 
The- significance  of  these  interactions  stresses  on  the  initiation  of  plastic 
flow  is  considered  and  it  is  shown  that  a  good  qualitative  agreement  has 
been  obtained  between  these  calculations  and  an  earlier  investigation  by 
Wojcik  and  Koss  (24).  The  effect  of  the  interface  phase  on  the  elastic  inter¬ 
action  stresses  has  been  considered  and  a  possible  role  of  interface  stresses 
on  elevated  temperature  creep  of  a  -  3  alloys  has  been  suggested. 


INTRODUCTION 


Increasing  interest  is  developing  in  the  role  of  colonies  of  Widman¬ 
statten  a  in  the  deformation  behaviour  of  a-0  alloys.  In  this  paper  we  ex¬ 
amine  the  conditions  leading  to  general  slip  across  a  colony  of  Widman¬ 
statten  a  surrounded  by  0. 

Slip  traversing  many  Widmanstatten  platelets  in  a  colony  in  various 
a-0  titanium  alloys  has  been  observed  by  several  investigators  including 
Wells  and  Sullivan(l),  Greenfield  and  Mangolin  (2),  Eylon  et.  al  (3),  Eylon 
and  Hall  (4),  Shechtman  and  Eylon  (5),  and  Eylon  and  Bania  (6).  Since  a 
precipitates  in  the  3  with  the  Burgers  orientation  relationship  (7) 

(000  H  { 110}  ^ 

<  1120  >a  H  <111  >3 

and  from  which  it  follows 

1  0“  10  }Q  //  {  lT  2 

<  11  20  //<  111  >3 

and 

{oil  l  }a  H  {  - 1  ”  )3 

<  11  20  >a  //  111  >3 

It  is  expected  that  slip  transfer  would  be  easy  for  these  Burger's  oriented 
slip  systems.  However,  this  orientation  relationship  does  not  explain 
what  is  the  effect  of  the  elastic  and  elastic-plastic  compatability  stresses 
at  the  a-0  interface  on  the  slip  in  a  and  3  and  across  the  colony.  It  will  be 
shown  later  that  long  range  elastic  interaction  stresses  exist  at  the  ar-3  in¬ 
terfaces  and  certainly  would  effect  the  slip  behavior  in  and  across  a  and  0. 
The  effect  of  the  interface  phase  (3.  Q,  10),  that  may  exist  in  some  of  the 


titanium  alloys,  on  the  slip  behavior  is  also  considered  and  will  be  dis¬ 
cussed. 

THEORETICAL  CONSIDERATIONS 

Hook  and  Hirth  (11)  have  shown  that  elastic  interaction  stresses  at  the 
grain  boundaries  of  bicrystals  give  rise  to  local  stresses  which  result  in  the 
operation  of  secondary  slip  systems  not  normally  expected  to  function  in 
the  individual  single  crystal  components.  Similar  behaviour  has  also  been 
reported  on  /3-brass  bicrystal  (12)  and  3-brass  bicrystal  and  tricrystal  (13). 
These  observations  can  be  looked  at  in  a  different  sense,  i.e.  the  inter¬ 
action  stresses  may  oppose  slip  in  some  slip  systems  and  aid  it  in  other 
slip  systems. 

If  one  assumes  a  cartesian  coordinate  system  xyz  such  that  x  and 
z  axis  lie  in  the  plane  of  the  two  phase  cr-3  interface  (Fig.  1).  Then  the 
continuity  of  strain  across  the  interface  is  related  by 


a  J3  a  _  .3  a  3 
e  *e  ,  =e. 


xx 


XX  zz 


zz'  ZX 


zx 


(1) 


The  other  strain  components  may  be  important  in  considering  the  effects  of 
end  constraints  but  do  not  affect  the  deformation  of  an  arbitrary  boundary 
(11).  It  is  possible  to  resolve  the  elastic  compatibility  strains  of  eqsTl] 
individually  onto  a  given  slip  system  and  add  them  to  determine  the  mag¬ 
nitude  of  the  resolved  elastic  strain  in  this  slip  system.  For  a  parallel 
slip  system  in  a  and  3  and  since  eqs  fl]  apply 

[2] 


a  3 

et  “  et 


where  efc  and  ^  are  the  total  resolved  shear  strains  from  eqs  [  1 1  of  a  and  3 
respectively. 

Let  us  consider  that  this  common  slip  system  is  a  primary  slip  system 
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of  <2  and  3.  Let  the  slip  plane  in  each  phase  be  the  yz  plane  and  the  slip 
directions  be  z,  Fig.  1.  The  orientation  of  this  slip  direction  is  not  the 
Burgers  orientation,  because  for  this  specific  orientation  the  slip  direction 
is  not  parallel  to  the  interface  as  Fig.  2  indicates. 

The  equations  relating  eqs[11  and[  2l  as  follows 

et=2e<*xx  coso^  cos +  2eazz  cos  03  cos  X. 3  +  ^zx  ^coa<^3'^ 

COS  X.J  TCOS  <X>  j  COS  \j) 

e^  =  2e^  cos  6.  cos  +  2e^  cos  o,  cos  \ 

t  xx  1  1  zz  3  i  1 

3 

+  e2X  (cos<t>3  cosX.^  +cos  6^  cos  \^)  ' 


where  and  <b  are  the  angles  between  the  slip  plane  normal  and  the  x  and 
z  axes,  respectively,  and  and  X.3  are  the  angles  between  the  slip  direction 
and  the  x  and  z  axis,  respectively. 

For  the  orientation  of  the  slip  plane  and  direction  given  earlier  and  to 
maintain  continuity 

a  _  a  3  _  3 


^  =  e  =  e;  =  e 


zx 


zx 


(4) 


Since  these  shear  strains  can  be  produced  by  a  single  shear  stress 


<rzx>  we  may  write  for  each  phase 


or  a  a  ,  ,  3  8  B 
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continuity  is  obtained  when  eq  [  5a  1  is  set  equal  to  eq  [  5b  1  .  If  e  /  e  ^  , 

then  interaction  stresses  will  develop  which  will  modify  the  magnitudes 

of  <r  and  <r  p  in  order  to  produce  equivalence  of  the  shear  strains. 
zx  zx 

Further,  if  e‘3  >  efl  then  3  will  tend  to  increase  e  and,  therefore,  ef- 

zx  zx  zx 

fectively  increases  the  shear  stress  on  the  a  slip  system. 


PROCEDURE 


A.  Interface  Plane  of  Alpha  and  Beta 

Martensite  in  Ti-Mn  alloys  has  been  shown  to  have  either  a  {  33-*}^ 
or  a  {  344  }^habit  plane  (14).  It  is  most  likely  that  martensite  and  Widman¬ 
statten  a  have  the  same  3  habit  plane  (15)  and  have  the  same  onientation  re¬ 
lationships  with  respect  to  the  parent  3.  Only  the  {  334  habit  plane,  how¬ 
ever,  has  been  reported  for  Widmanstatten  a  (9,  16,  17),  and  consequently 
{  334  }gis  taken  to  be  the  plane  of  3  parallel  to  the  interface.  If  the 
Burgers  orientation  relation  is  arranged  so  that 

(0001)^  //  (110)^ 

[  27  To  }a  //  [iTTJ3 

then  the  plane  of  a  parallel  to  (  334)^  is  (  13  580  )  as  is  shown  in  Fig.  2. 

If,  however,  the  [12  10  1  ^  direction  is  made  parallel  to  the  [11  1 »  which 
also  lies  in  (110) 0  ,  then  the  plane  of  a  parallel  to  1 3 3 4 ) ^  is  (5140).  The 
existence  of  the  {5140  }  parallel  to  {334}o  was  originally  reported  by 
Albert  (18)  and  in  recent  investigations  (9,  17).  The  pole  of  the  (5140)  plane 
is  shown  in  Fig.  2  and  it  is  11°  away  from  (13580). 

Most  of  the  calculation  reported  here  have  been  made  for  the  (  13  580  ) 
interface  plane.  It  will  be  shown  that  at  the  point  when  there  is  maximum 
interaction  between  alpha  and  beta  the  difference  between  the  elastic 
strains  on  the  slip  systems  when  calculated  for  the  two  interface  planes  of 
alpha  is  about  5%. 

B.  Elastic  Constants 

The  five  elastic  compliances  for  single  crystals  of  a  (19)  are  s^  = 
0.997,  0.472,  s^  =  -  0.193,  s^=0.  688,  s^4  =2.14  in  units  of  10 

cm4 /dyne  .  Fisher  and  Dever  (20)  have  measured  the  elastic  stiffnesses  of 


3-Ti  -  10  wt.  7i  cr  crystals  at  25'' C  as  =  1.  331,  c,0  =  0.  951  and  = 

12  o 

0.427.  In  units  of  10  dynes/cm".  When  converted  to  elastic  compliaxt- 

-12 

ces,  the  values  are  =  1.  857,  s^=  "0*  7~4  and  =  2.  342  in  units  of  10 

cm"/ dyne.  These  compliances  refer  to  the  principal  crysrallographic 

ci  a  a  3  3  3 

axes,  indicated  by  squares  in  Fig.  2,  xq,  y^,  and  x^  ,  yQ  and  z  ^  . 

In  the  calculations  made  a  series  of  axes,  at  various  positions  with 
respect  to  the  interface  were  assumed.  It  was  necessary  to  convert  the 
elastic  compliances  to  each  of  these  axes.  The  necessary  transformation 
equations  can  be  derived  according  to  procedures  outlined  by  Zener  (21) 
or  Nye  (22). 

If  one  assumes  that  only  a  normal  stress  is  applied  on  a  crystal  in  a 
given  direction  and  if  this  direction  is  set  to  be  x'(y'  and  z'  is  obtained  so 
that  x'y'z'  will  form  a  right  handed  coordinate  system)  then  to  define  the 
state  of  strain  in  the  crystal  it  is  necessary  only  to  calculate  the  six 
elastic  constants  namely  s'^,  s'^,  s  41  s ' 51  and  to  find  the  six 

strains  as  shown  below; 


»  ,  ,  =  s',.  o-  ,  .,e,  , 
xx  11  x  x  y  y' 


s  .,.  or  .  .  ,  e  .  , 
21  xx  z  z 


s  '  ,  <r  .  , 
31  xx 


'y'z' 


=  s '  <t  ,  ,  e  ,  ,  =  s',,  <r  ,  ,  and  e  .  ,  =s',.  cr  ,  . 

-tl  xx  z  x  51  xx  x  y  61  x  x 


(6) 


where  <r  ,  ,  is  the  applied  stress  and  is  considered  constant.  The  six 
elastic  constants  were  calculated  for  both  phases  for  a  total  54  different 
directions  ,  Fig.  2. 

C.  Elastic  Strains  Resolved  into  Interface  and  Slip  Systems 

Once  the  six  elastic  compliances  for  each  stress  direction  for  a  and  3 
have  been  calculated,  the  strains  could  be  obtained  from  eq.  [  61  .  The  next 
step  was  to  transform  these  strains  into  a  new  set  of  axes  xyz  such  that  xz 


■  '.jr'TrrvnQi*. 


is  the  interface  plane  as  indicated  in  Fi  gs.  1  and  2.  Since  the  transformation 
laws  for  such  a  tranformation  of  strains  is  the  same  for  both  the  a  and  3 
phases  one  can  avoid  transforming  the  strains  of  all  the  54  stress  directions 
into  the  new  set  of  axis  xyz  by  examining  the  strains  to  obtain  represent¬ 
ative  results.  The  elastic  strains  of  the  a  and  3  phases  for  each  stress 
direction  of  the  54  stress  directions  were  compared  and  the  stress 
directions  for  the  transformation  of  strains  were  chosen  so  that  in  some 
stress  directions  maximum  and  in  other  directions  minimum  interactions 
was  likely  to  occur. 

The  chosen  stress  directions  namely  x x'^  «tc.  are  circled  as 
shown  in  Fig.  2.  The  strains  in  x'j  y'j  z  j  ,  in  x'j,  y  z'5  etc.  sets  of 
axes  were  then  transformed  into  the  new  set  of  axis  x,  y.  z.  Figs.  1  and  2. 

As  indicated  earlier,  eq.  [11,  the  strains  which  need  to  be  considered  for 
compability  are  e  ,  e  and  e  ,  and  these  elastic  strains  for  both  the 

XX  Z  Z  Z  X 

n  and  3  phases  have  in  turn,  been  resolved  into  12  oe  slip  systems,  3  prisom 
3  basal  and  6  pyramidal.  It  was  assumed  chat  the  a  slip  systems  had 
parallel  counter  parts  in  the  3  phase,  although  these  planes  and  directions 
would  not  necessarily  act  as  slip  systems  in  the  3,  and  elastic  strain  of  3 
were  resolved  onto  these  'slip  '  systems.  The  total  resolved  compatibility 
shear  strains  on  any  given  slip  system  were  obtained  from  eq.  [  3  ]. 

D.  Maximum  Interaction  Stress 

An  attempt  was  made  to  determine  approximately  the  interaction  stress 

corresponding  to  the  stress  direction,  x'  for  which  the  difference  between 

a  3 

et  and  e^for  the  slip  system  of  a  and  3  was  found  to  be  maximum.  For  cr 
this  slip  system  is  designate  E,  (0001)  [1210]  ,  Table  I.  The  corresponding 
"slip"  system  of  beta  is  not  of  specific  interest.  The  shear  strains  e^  for 


slip  system  E  and  the  corresponding  e{  were  determined.  The  corresponding 
values  of  the  elastic  compliances  s  ’  and  sH’  were  calculated  as  2.141 
and  3.415  in  units  of  lO'1^  cm2/dyne,  respectively.  Eastic  shear  stress-shear 
strains  curves  for  E  slip  systems  are  shown  in  Fig.  4.  From  these  two 
curves  it  is  possible  to  make  a  first  order  calculation  of  the  interaction 
stresses,  as  well  be  discussed  in  the  results,  by  assuming  that  a  stress  of 
a  specific  magnitude  is  applied  along  the  x'g^  direction. 

E.  Interface  Phase 

Although  the  volume  fraction  of  interface  phase  is  quite  small  and  its 
thickness  is  usually  less  than  0.4  -m  (IT),  some  consideration  has  been 
given  to  its  possible  role  on  producing  interaction  stresses  because  of  the 
possibility  that  it  could  have  a  high  yield  strength  and  therefore  withstand 
high  total  elastic  stresses.  When  the  interface  phase  is  alpha  (rather  than 
the  transitional  fee  phase  (17)  ),  it  has  been  reported  to  be  either  in  {  10  1  1} 

(10)  or  { 10  1*2}  (8,  9)  twin  relationship  with  respect  to  the  primary  a.  The 
{10  1  l}twins  have  a  Burgers  orientation  relationship  with  3( 23)  and  hence 
these  specific  {10  11}  twins  have  been  considered  in  this  anaylsis.  It  was 
not  possible  to  ascertain  which  specific  {10  12}  twins  were  present  and  ac¬ 
cordingly  all  six  possibilities  were  considered.  Fig.  3.  Interactions  have 


been  determined  only  for  the  stress  direction  x'24.  The  relative  posi¬ 
tions  of  the  xaQ,  y*  and  z°^  change  with  respect  to  the  x1  y'  4>  z' ^ 

coordinate  system.  The  x  y  ^ are  the  new  positions  of  the  xQ  ,  y° , 

OOO  OO 


q  Ct  Q  Ct 

zq  positions  for  the  twin  G(0  111)  .Once  the  new  positions  of  x^,  yQ.  z^ 


have  been  determined  then  the  procedure,  outlined  previously  for  the  deter¬ 
mination  of  interactions,  was  followed. 


RESULTS 


A.  Total  Resolved  Compatibility  Shear  Strains  » 

The  elastic  compliances,  and  hence  the  strains,  have  been  calculated 
for  the  54  stress  directions.  For  the  twelve  directions  of  interest  the 
strains  have  been  transformed  to  the  x,  y,  z  axes  related  to  the  interface, 

Figs.  1  and  2,  and  then  the  e  ,  e  and  e  strains  were  further  resolved 
on  to  the  twelve  slip  systems,  previously  given,  for  both  the  a  and  3  phases, 
according  to  eq.  F  3  1 .  The  results  are  given  in  Tables  la  and  lb.  The 
Schmid  factors  are  given  in  Table  2.  The  Schmid  factors  have  positive  and 
negative  signs  because  positive  directions  have  been  assigned  to  the  slip 
directions  irrespective  of  the  applied  stress.  The  last  column  of  Table  1, 
contains  comments  indicating  whether  the  interaction  strains  aid  or  oppose 
the  applied  strains  induced  directly  from  the  applied  stress.  The  interact¬ 
ion  strains  aid  or  add  to  the  applied  strains  when  the  difference  between  the 
strains  in  3  and  a,  Ae  =  e^  -  e^,  has  the  same  sign  as  the  applied  strain. 

The  interaction  strains  oppose  the  applied  strains  when  this  difference  has 
a  sign  opposite  to  the  applied  strain. 

In  Tables  la.lb,  and  3  the  x'*4  direction  is  for  the  case  where  the  in¬ 
terface  planes  are  (T  140)o  and  (T  34)^.  For  this  pair  of  interface  planes 
the  Burgers  orientation  is 

(0301  )aH  UlO)^ 

[i2ioia  //  [ n  i] £ 

It  is  to  be  noted  that  the  x'*4  directions  is  similar  to  the  x'24  direction,  i.  e. 
both  directions  lie  in  their  respective  interface  planes  and  are  both  50°  away 
from  [  0001]a  and  [110  . 
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B.  Total  Resolved  Compatibility  Shear  Strains  for  the  Interface  Phase 

Table  3  gives  the  total  resolved,  compatibility-  shear  strains  resolved 
onto  the  twelve  slip  systems  of  a,  designated  A,B  etc.  The  calculations 
were  carried  out  in  the  same  manner  as  for  Table  1,  for  the  direction  x1 

■C,  Estimate  of  Compatibility  Shear  Stresses 

When  a  given  stress  is  applied  in  a  given  direction  in  a  and  3*  the 
strains  when  resolved  into  the  interface  will  generally  not  be  the  same.  Ac¬ 
cordingly  some  adjustment  must  be  made  in  the  interface  strains  of  a  and 
3  with  the  result  that  interaction  stresses  will  arise,  as  pointed  out  earlier. 

When  the  adjustments  have  been  made,  all  the  compatibility  strains 
will  be  equivalent  and  when  resolved  onto  the  parallel  slip  systems  of  or 
and  3  the  elastic  strains  on  these  slip  systems  must  also  be  equal. 

Let  us  assume  that  a  constant  stress,  say  103.4  MPa  (15  ksi  or  10.3.4 

x  107  dynes/cm^)  is  applied  to  each  phase  along  the  x' direction*  This 

o  E  8  E 

stress  will  produce  the  strains  efc  ’  and  e^  ’  ,  see  Table  1.  If  one  assumes 

v  E  3  E 

that  equivalent  shear  stresses  < r  ’  and  <r  '  have  been  applied  to  pro- 
duce  these  strains,  then  the  e^’  and  e£*  positions  will  be  located  as  shown 
in  Fig.  4.  It  is  assumed  that  there  are  equal  volume  fractions  of  each  phase. 
Sine  the  shear  stresses  in  each  phases  are  parallel  and  the  cross  sectional 
areas  are  equal  a  decrease  in  stress  in  one  component  must  be  matched  by 
an  increase  in  stress  in  the  other.  Thus  we  may  write 


A  or,  E  A  3,  E 
A  <r  =  A  <r 


(7) 


A<r  a'E  =  (e"’E 


e*’E  )-i- 

t  VlE 


(8) 


A  tr 


3,  E 


3,  E 


.3,  E 


1 


(9) 
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Where  Air E  and  Atr^'E  are  the  changes  in  equivalent  shear  stress  along 

,  ,  a,  E  ,  8.  E 

the  E  slip  system  for  a  and  3,  respectively,  e£  ^  and  are 

the  final,  equal  shear  strains  on  the  E  slip  systems  after  compati¬ 
bility  has  been  established. 

From  eqr  s  [71,(8]  and  [  9  1  it  follows  that 


<*.  E  s3,  E  +  J,  E  a.  E 
a,  E  3,  E  et  •  S  +  et  •  3 
6t,f  =Ct.f  =  t«.'E  * 


(10) 


The  strains  e0f’E  and  ef*  E  are  shown  on  their  respective  curves 
t,  f  t,  t  /  2 

.  .  ,  „  lrt7  dynes/ cm 

in  Fig.  4.  The  change  in  stress  is  14  MPa  (2.0ksi,  14x10 

which  is  30%  of  the  resolved  shear  stress.  This  increment  of  stress  in  a 

adds  to  the  applied  stress  and  thus  increases  the  total  stress  on  the  basal 

slip  system  E. 


DISCUSSION 

A.  Effect  of  Stress  Direction  on  Elastic  Interactions 

The  interaction  stresses  are  only  a  fraction  of  the  shear  stresses 

produced  by  the  applied  stress.  Consequently,  these  stresses  will 

be  significant  in  initiating  slip  only  for  those  slip  systems  for  which  the 

Schmid  factors  are  highest.  The  interactions  are  highest  for  the  stress 
* 

directions  x'5g  ,x'24  and  x'^,  Table  la  and  lb,  Fig.  2.  For  example,  for 
stress  direction  x'24  the  total  resolved  compatibility  elastic  strain  on  the 
E  slip  system  of  a  and  3  are 

<’E  =  0.978  o- 
1  x  24 

e/3,  E  =1.718  <r  ,  . 
t  x  24 


A.  ..f'E 


ef'  E  =  0.740  e  ... 
t  x  24 
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The  large  value  of  de  indicates  that  the  elastic  strins  in  0  will  strong 


ly  assist  slip  in  a,  since  the  signs  of  A  e  and  the  Schmid  factor  are  the 
same.  The  largest  value  of  A  e  is  Ae  =0.90?  o-^.g  for  the  E  slip  system 
and  here  £  will  also  assist  slip  in  a. 

Interactions  can  also  be  quite  high  for  x'gg  and  x'^,  but  *or  these 
slip  systems  the  Schmid  factors  are  low. 

The  increases  to  the  resolved  shear  stress  of  a  by  the  elastic  inter¬ 
action  was  estimated  to  be  30%  for  x'0^,  Fig.  4  on  the  E  slip  system,  where 
the  maximum  interactions  are  found.  Examination  of  Tables  la  and  lb  in¬ 
dicates  that  the  interactions  on  the  E  slip  system  are  somewhat  higher  for 
the  x'*^  stress  direction,  which  lies  on  the  (5  140  )  interface  plane  paral¬ 
lel  to  (334)g.  A  calculations  analogous  to  that  for  x1  ^  was  also  carried  out 

for  the  x'gg  directions,  and  the  increment  to  the  resolved  shear  stress  of 

* 

a  was  35%  »  the  maximum  increment.  The  x'  ,  x  „  .  and  x'_0  axes  all 

44  44  DO 

lie  near  the  (13580)^  interface  plane.and  lie  between  40  and  50°  to  the  pole 

of  the  basal  plane.  The  x'  and  x',.,,  axis  also  lie  within  40-50°  of  the  basal 

4  i  d  4 

plane  pole,  but  along  a  different  plane.  Fig.  2. 

Within  a  range  of  10°  to  the  basal  plane  pole,  stress  axes,  x  ^g.  x'  .^q> 
x'-g.  Fig.  2,  the  interactions  are  small,  Tables  la  and  lb.  When  the  stress 
axis  lie  in  the  basal  plan,  x  x'g,  ,  x  x'^q.  Fig.  2,  there  is  no  inter¬ 

action  for  the  basal  slip  system.  For  these  stress  axes  some  of  the 
Schmid  factors  for  the  prism  and  pyramidal  slip  are  high.  In  these  cases 
the  interaction  stresses  weakly  oppose  the  applied  stress.  When  compared 
to  other  slip  systems  interactions  are  strongest  for  the  basal  slip  systems 
and  thus  slip  on  the  basal  slip  systems  can  be  most  strongly  aided  by  these 
stresses.  Elastic  interactions,  although  weak,  can  also  assist  pyramidal 
slip  ,  see  slip  system  U for  stress  axis  x!.?^,  Table  lb.  However,  for 
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prism  slip,  elastic  interactions  only  weakly  oppose  the  strains  resulting 
from  applied  stress,  see  slip  system  C  for  stress  axis  x'j,  Table  la. 


1,  Role  of  Interaction  Stresses 

In  the  following  discussion  it  is  assumed,  for  lack  of  data,  that  the 
elastic  constants  of  pure  a  (19)  and  3  containing  10  wt  %  cr  (20)  apply  for  a 
and  3  phases  of  Ti-  3  Al-  IMo  -IV,  results  from  which  alloy  (24)  are  to  be 
discussed  here.  The  elastic  interactions,  which  have  been  presented 

.  Tf 

earlier,  take  place  at  a  single  a-B  interface.  In  a  colony  of  Widmanstatten 
a  and  3  ,  the  or  ientations  of  both  a  and  3  are  constant.  Thus,  when  sig¬ 
nificant  interaction  stresses  are  present,  they  occur  at  each  interface  and, 
thus,  a  long  range  elastic  interaction  stress  is  created. 

Wojcik  and  Koss  (24)  have  determined  the  resolved  shear  stress  on 
the  observed  slip  plane  at  the  0.2%  yield  stress  for  colonies  of  Widman¬ 
statten  or  and  3  in  a  Ti-8  Al-lMo-IV  alloy.  The  orientations  of  the  stress 
axis  in  a  standard  stereographic  triangle  are  presented  in  Fig.  5  and  the 
results  are  given  in  Table  4.  It  can  be  seen  that  there  is  a  range  of  resolved 
shear  stresses  from  228  MPa  for  stress  axis  15  to  463  MPa  for  stress  axis 
17. 

To  assess  the  role  of  the  interface  interactions  in  affecting  these 
stresses  the  orientation  of  the  stress  axis  with  respect  to  the  a-3  interface 
must  be  considered.  For  stress  axis  15  the  pole  of  the  3  interface  plane  is 
95°  from  the  stress  axis  and  for  stress  axis  17  the  3  interface  pole  is  77° 
away  (24).  From  this  information,  and  that  of  Fig.  5,  axis  15  is  located  at 
x'^g  and  axis  17  is  located  at  x'-g  in  Fig.  2. 

Since  the  applied  stresses  are  known  for  these  stress  axes,  Table  4, 
calculations  similar  to  those  made  for  x'94  were  carried  out  to  ascertain 
the  approximate  interaction  stresses.  For  stress  axis  x'^g  the  elastic 


in  Colonies  of  Widmanstatten  a  and 


strain  in  3  is  higher  than  the  strain  in  or,  and  if  the  volume  fractions  of  a 

and  3  were  equal,  90  MPa  would  be  added  to  the  nominal  stress  for  slip  in 
* 

Table  4  for  slip  system  E  with  the  highest  Schmid  factor.  For  stress 
axis  x',g,  3  opposes  a  deformation  and  again  for  equal  volume  fractions  of 
a  and  3,  the  decrease  in  nominal  stress  would  be  29  MPa  for  slip  system 
F,  for  which  the  Schmid  factor  is  highest.  It  should  be  noted  that  when 
the  interaction  stress  aids  slip  the  stress  required  from  the  applied  stress 
is  reduced,  and  when  the  interaction  stress  opposes  slip  a  greater  applied 
stress  is  required. 


It  is  clear  that  when  the  elastic  interactions  for  stress  axis  x'gg  and 
x'^g  are  considered  the  differences  between  resolved  shear  stresses  for 
slip  for  axis  15  and  17  would  be  reduced.  If  the  maximum  interactions  took 


place,  then  the  stresses  for  slip  would  be 


axis  15  (x'5g)  : 

363 

MPa 

axis  17  (x'59)  : 

485 

MPa 

The  difference  in  resolved  shear  stress,  when  elastic  stresses  are  con¬ 
sidered,  would  be  about  34%  compared  to  axis  15  stress  and  are  88. 4%  when 
the  elastic  stresses  are  not  considered.  One  may  reasonably  question 
whether  the  maximum  interaction  is  possible  since  the  volume  per  cent  of 
3  is  considerably  less  than  50,  nevertheless,  the  trend  must  exist. 

It  also  appears,  reasonable  to  believe  that,  when  the  0.2%  yield  stress 

is  reached,  slip  in  3  as  well  as  a  is  taking  place.  If  slip  in  a  is  initiated 
* 

Given  the  orientations  shown  in  Fig.  5,  calculations  as  to  the  resolved  shear 
stress  for  axis  15  (x'gg)  and  17  (x'g„)  are  somewhat  different  from  those 
given  in  Table  4.  They  are  respectively  273  and  514  MPa  instead  of  228  and 
463  MPa. 
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before  slip  in  3,  then  this  slip  in  a,  which  will  pile  up  dislocations  at  the 
interface,  will  help  to  initiate  slip  in  3.  This  assistance  would  contribute 
to  the  lower  resolved  shear  stress  for  axis  15.  The  stress  axis  *'53 
(axis  15)  and  x'-^  (axis!7  )  have  been  referred  to  the  (13580)^  interface 
plane.  There  is  little  difference  in  the  results,  if  they  are  oriented  with 
respect  to  the  (Tl40)ff  interface  plane.  It  is  of  interest  to  note  that  the 
lowest  observed  resolved  shear  stress,  axis  15,  occurs  when  the  assistance 

to  slip  in  a  by  3  is  highest. 

Table  4,  shows  that  other  slip  systems  may  act  in  additions  to  basal 
slip.  For  appropriate  oreintation,  the  Schmid  factor,  and  elastic  inter¬ 
actions  may  favor  slip  in  these  systems.  It  has  been  pointed  out  that  the 
<r-3  interfaces  are  not  always  planar  (24)  and  that  this  may  contribute  to 
the  differences  in  stress  noted  in  Table  4.  Further,  other  factors  now 
being  considered  (24),  may  reduce  still  more  the  differences  in  stress 
which  remains  unaccounted  tor  by  elastic  interactions. 

2.  Role  of  the  Interface  Phase 

For  the  stress  direction  x'24  calculations  were  carried  out  for  8  dif¬ 
ferent  kinds  of  twin  interface  orientations,  Fig.  3.  The  calculated  elastic 
strains  are  given  in  Table  3.  To  evaluate  the  effect  of  the  total  resolved 
compatibility  strains  of  the  twin,  e£,  one  must  compare  these  strains  with 
the  corresponding  strains  in  a  for  the  x' 24  ax^-s>  Tables  la  and  lb  for  all 
12  slip  systems  of  or  and  the  corresponding  12  parallel  slip  systems  of  the 
twin.  The  comparsion  will  be  made,  for  the  moment,  independent  of  the 
relative  sizes  of  the  two  regions. 

From  Tables  la,  lb  and  3,  it  can  be  seen  that  the  differences  between 
e^  and  ek  ,  in  most  instances,  are  very  small.  For  example,  compare  the 

e^’  for  x'24,  Table  la,  with  e  ’  for  the  J  twinning  plane,  .Table  3: 
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e®’  =  Q.9"8  ir  and  e|’  ^  =0.953  <r  x'24  ,  a  difference  of 

0.025  <r  x'24. 

When  there  are  significant  differences,  then  the  strains  of  the  twins 

are  in  between  the  strains  of  the  a  and  3  phases.  For  e®’  ^  ,  e^’  ^  and 
t  E 

et’  for  the  N  twinning  plane  of  Tabie  3,  the  values  are  as  follows. 

E  •  878  *  x  24 

et’  E  r  1*  299  <r  x'24 

%Bl  E  =  1.  7>8  x  *'24 


In  these  circumstances  the  twin  acts  in  the  same  direction  as  3,  assisting 
deformation  in  a. 

Only  for  the  G  twin  does  the  twin  oppose  deformation  of  a  and  here 
the  magnitude  of  the  opposition  is  small.  Compare  the  following  et  values. 


«“’E  .  0.978  <r*'24  | 

*'’E  =  0.831.  x'24  I 


0.147  <r  x’24 


e®,E  =  1.718  O'  x’24 

However,  the  twin  strongly  opposes  the  strain  in  3  and  the  net  effect  is  that 
the  assistance  to  slip  in  or  by  3  is  reduced.  Thus  interactions  similar  to 
those  found  at  <v/3  interfaces  can  be  found  at  twin/3  and  at  a/twin  interfaces. 

The  thickness  of  the  twins  can  extend  up  to  about  0.4jum(17)  which 
is  much  smaller  than  the  thickness  of  a  and  3.  Thus  the  role  of  twins  in 
assisting  or  resisting  deformation  of  o  by  3  is  likely  to  be  quite  limited. 
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Crystallographically  the  twins  at  the  interface  can  readily  permit 
transfer  of  slip  from  or  to  twins  to  B,  because  of  the  existence  of  the  follow¬ 
ing  relationships,  Fig,  3, 


r 


{l07l}frwifl  G,  H 


(oooi>_  a  {  )//  {no}  B 

or  H 

^10  1  twin  I,  J,  K,  L,  M,  N(within  5°) 


<  11  2  0  >  //  <  11  2  0  > 


a 


twin  G,  H,  I,  J, 


//  <  111  >, 


for  twin  planes  M  and  N  transfer  of  slip  to  B  should  be  some  what  more 
difficult,  because  the  <  11  20  >  direction  of  both  a  and  twin  are  about  11°  from 
<11  1>3.  For  the  K  and  L  twin  planes  transfer  of  slip  from  the  twin  to  B 
c-n  not  easily  be  achieved  because  of  the  large  angle  between  the  slip  direc¬ 
tions,  Fig,  3. 

For  the  E  slip  system, (0001)  [  1  2  10  ]  ,  for  which  the  elastic  interactions 
are  greatest, Fig.  3  indicates  that  [1  210  ]  is  aboutli""  away  from  [1  1  1]^ 
for  M  and  N  {  10  12}  twins, and  thus,  transfer  of  slip  would  be  intermediate 
in  difficulty.  However,  if,  in  setting  up  the  Burgers  orientation  relation¬ 
ship,  these  two  directions  had  been  made  parallel, the  a  habit  plane  would 
have  become  (5  140)  and  the  magnitude  of  the  elastic  interactions  would  have 
been  about  the  same  as  for  the  (  1  3  580  )  interface,  as  previous  calculations 
have  shown.  In  this  case,  slip  tranfer  would  have  been  quite  readily  achieved 
for  both  (10  11  }  and  (10  12  }  twins.  Thus,  depending  on  which  twinning 


system  occurs,  it  is  quite  possible  for  the  crystallography  to  permit  ready 
transfer  of  slip. 

For  those  instances  where  transfer  of  slip  would  be  easy,  only  the 

dislocations,  and  defect  structure  of  the  twin,  developed  during  its  for¬ 
mation  would  serve  as  an  obstacle  to  slip  transfer. 

C.  Elevated  Temperature  Effects  of  Interface  Stresses 

In  the  preceding  discussion  on  interaction  stresses  it  has  been  pointed 
out  that,  at  room  temperature,  &  can  assist  slip  in  a,  as  a  result  of  inter¬ 
action  stresses.  At  elevated  temperatures,  it  has  been  reported  that 
WidmanstUtten  structures  have  higher  creep  resistance  than  equiated  struc¬ 
tures.  It  has  long  been  recognized  that  the  a  phase  has  higher  creep  re¬ 
sistance  than  the  j3  phase.  A  change  in  morphology  does  not  change  the  com¬ 
position  of  the  phases.  Consequently,  one  may  look  to  the  difference  on 
creep  behavior  of  the  two  morphologies  in  terms  of  the  structure. 

One  simple  way  of  looking  at  the  situation  is  to  consider  that  the 
existence  of  the  Widmanstatten  structure  reduces  the  path  lengths  of  slip 
on  the  j3  and,  as  a  reslult,  increases  the  creep  resistance. 

Another  point  of  view  is  that  elastic  interaction  stresses  may  contri¬ 
bute  to  creep  resistance.  Although  ,  3  may  assists  deformation  of  a,  or  at  the 
same  time  would  impede  deformation  of  3.  At  elevated  temperatures, 
when  3  has  les3  resistance  to  creep,  the  resistance  to  deformation  in  3. 
provided  by  a,  would  increase  creep  resistance.  Further,  since  the  inter¬ 
action  stresses  occur  over  an  entire  colony,  the  whole  region  would  be  ex¬ 
pected  to  be  strengthened  by  the  interaction  stresses.  In  addition  to  the 

inherent  behavior  of  each  phase  and  the  effect  of  interaction  stresses,  one 
would  also  expect  that  volume  f  vact  ion  of  phases  in  the  Widmanstatten  co¬ 
lony  structures  would  also  be  significant. 
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SUMMARY 


1.  The  six  elastic  strains,  as sociated  with  stress  along  a  single  nor¬ 
mal  axis,  have  been  calculated  for  54  different  stress  axes.  These  strains 
were  then  resolved  along  the  a-3  interface  for  twelve  representative  direc¬ 
tions.  Only  those  strains  associated  with  compatibility  were  resolved  onto 
twelve  slip  systems  of  a  and  a  corresponding  set  of  parallel  slip  systems 
in  3.  The  slip  systems,  assumed  for  3,  were  not  necessarily  actual  slip 
systems.  From  a  comparison  of  these  resolved  elastic  strains  it  was  pos¬ 
sible  to  determine  whether  the  elastic  strain  in  3  would  add  to  the  elastic 
strain  in  a  and,  thereby,  assist  in  the  initiation  of  slip  in  a. 

2.  Calculations  were  carried  out  for  the  Burgers  orientation,  in  which 

moi)a  //  {no 

<  n  20  >  //  <  iu>  - 

a  4 

and  for  the  assumption  that  the  interface  plane  of  a  parallel  to  {3  34  was 
{  13  530  }a  .  The  alternate  interface  plane  of  a,  {5  140  }Q  ,  reported  in  the 
literature  occurs  when  the  second  of  the  <  111  >g  directions  in  {110}^  is 
chosen.  Calculations  shows  that,  if  {5  140}^,  is  the  interface  plane,  there 
is  very  little  difference  in  the  magnitude  and  type  of  interaction  stress. 

3.  Maximum  elastic  interactions  between  a  and  3  were  found  when  the 
stress  axes  were  in  or  close  to  the  interface  plane.  An  estimate  of  the  mag¬ 
nitude  of  these  stresses,  based  on  compatibility  strains  alone,  indicated 
that  the  elastic  interaction  stresses  could  be  30-357o  of  the  resolved  shear 
stress. 

4.  It  was  found  that  3  provided  maximum  interaction  assistance  to  strain 
in  a  along  the  basal  slip  system  (0001)  [  1210  1.  Interaction  assistance  was 
found  to  be  smaller  for  prism  and  pyramidal  slip. 
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5.  Calculations  for  elastic  interactions  in  colonies  of  Widmanstatten  a 
and  3  were  carried  out  for  the  alloy,  Ti-8Al-lM0-IV,  for  which  experi¬ 
mental  data  were  supplied  by  Wojcik  and  Koss  (24).  It  was  found  that  the 
highest  interaction  assistance  from  .3  occurred  for  the  basal  slip  systems 
for  which  the  experimental  resolved  shear  stress  was  least. 

6.  It  was  proposed  that  the  elastic  stresses  at  a-/3  interfaces  of  Widman- 

ii 

statten  structures,  stresses  which  oppose  slip  in  p,  contribute  to  the  greater 
creep  resistance  of  Widmanstatten  over  equiaxed  a- 3  structures  at  elevated 
temperatures . 
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Fig.  1. 
Fig.  2. 


Fig.  3. 


FIGURES 

Geometry  of  the  Alpha-Beta  interface  and  coordinate  axes. 


A  stereographic  projection  showing  the  interface  plane, 

(13580)  ,  (334)^  ,  Burgers  orientation  relationships,  principal 

e  °  e  .  ..  .  ar  or  ar  ,  p  ft  ft 

reference  axes  for  elastic  constants,  x  ,  y  ,  z  and  x  ,  y  ,  z  , 

o  7o  o  o  7o  o 


and  3tress  directions,  xi 


A  stereographic  projection  showing  the  poles  and  planes  of  the 

{lOll  }  twins,  G.  H  and  the  {10T2  }  twins,  I,  J,  K,  L,  M,  and 
N;  and  the  crystallographic  relationships  between  or  and  the 
twinned  cr. 


Fig.  4.  The  determination  of  the  magnitude  of  the  interaction  stress  for 
the  stress  direction  x^on  the  slip  system  E. 

Fig.  5(24).  The  orientation  of  the  stress  axes  in  a  standard  HCP  stereo- 
graphic  triangle  for  a  in  Widmanstatten  colonies. 
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Resolved 


•mi»: 


TABLE  4  (24) 

Compressive  Yield  Stress  end  Resolved 
Sheer  Stress  on  the  Observed  Slip  Planes. 
Sample  quenched  from  760 °c  unless  other¬ 
wise  noted. 


Sample 

Number 

0.2%  Yield 
Stress 

Resolved  Shear  Stress 
on  Observed  Slip  Plane 

6 

1160  MPa 

414  MPa 

7 

634 

312 

8 

2190 

245 

9* 

612 

281 

10+ 

731 

348 

11 

876 

341 

12*+ 

724 

344 

13+ 

1610 

394 

15 

554  * 

228 

17 

1890 

463 

19+ 

1480 

454 

20 

605 

235 

21 

590 

262 

22* 

636 

251 

♦These  samples 

did  not  slip  on  the 

(0001)  plane. 

tHeat  treated  at  925  9C  and  quenched. 
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FIGURE  3 
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PART -II;  FINITE  ELEMENT  METHOD  (FEM)  CALCULATIONS  OF  STRESS- 
STRAIN  BEHAVIOR  OF  ALPHA- BETA  Ti-Mn  ALLOYS 

ABSTRACT 

By  use  of  a  NASTRAN  (18)  Computer  Program,  the  Finite  Ele¬ 
ment  Method  (FEM)  has  been  employed  to  calculate  the  effect  of 
particle  size,  matrix  and  volume  fraction  on  the  stress-strain  relations 
of  or— p  titanium  alloys.  It  was  found  that  for  a  given  volume  frac¬ 
tion,  the  calculated  stress-strain  curve  was  higher  for  a  finer  particle 
size  than  for  a  coarse  particle  size  within  the  range  of  the  strains 
considered,  and  this  behavior  was  seen  for  all  the  different  volume 
fraction  alloys  considered.  For  a  50:50  vol  %  ct-p  alloy,  the  stress- 
strain  curve  with  p,  the'  stronger  phase,  as  the  matrix  was  higher 
than  that  with  a,  the  softer  phase,  as  the  matrix.  The  calculated 
stress-strain  curves  for  four  different  vol.  %  a  alloys  were  compared 
with  their  corresponding  experimental  curves  and  in  general  good 
agreement  was  found.  Whenever  there  were  discrepancies,  they  were 
discussed  by  comparing  the  morphology  of  the  mesh  used  in  the 
calculations  with  the  morphology  of  the  actual  materials. 
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Two  ductile  phase  structures  occur  in  many  technologically  im¬ 
portant  alloys.  Increasing  interest  is  developing  in  calculating  the 
stress-strain  behavior  of  the  two-phase  structures  from  the  stress- 
strain  curves  of  the  component  phases.  In  the  present  calculations 
FEM  has  been  employed  to  calculate  the  eifecr  of  particle  si2e,  matrix 
and  volume  fraction  of  phases  on  the  stres:  -strain  behavkr  of  a  -0 
Ti  alloys.  In  Part  I,  stress -strain  relations  will  be  shown  and  com¬ 
pared  with  the  experimental  curves  for  4  differenr  volume  fractions. 
In  Part  II,  stress-strain  distributions  will  be  shown  and  strain  distri¬ 
butions  will  be  compared  with  experimentally  determined  strain  distri¬ 
butions. 

Two  simple  hypotheses  have  been  suggested  by  Dorn  ana 
Starr  (1)  to  calculate  the  stress -strain  curve  of  a  two  ductile  phase 
mixture  from  the  stress-strain  curve  of  each  phase  and  volume  frac¬ 
tions:  one  is  based  on  constant  strain  in  each  phase  and  the  other  is 
based  on  constant  stress  in  each  phase.  Usually,  the  experimental 
curves  of  the  two  phase  mixtures  will  lie  somewhere  between  the  two 
curves  which  are  calculated  according  to  the  above  two  hypotheses,  a 
fact  that  Dorn  and  Starr  anticipated. 

Davies  (2,3,4)  studied  the  deformation  behavior  of  dual  phase, 
ferrite-martensite,  steels  and  have  found  that  the  0.2%  YS  linearly 
varied  with  the  volume  fraction  of  the  martensite,  and  followed  a  law 
of  mixture  rule.  Interestingly,  he  (2,3,4)  also  found  that  the  tensile 
strength  and  ductility  of  the  dual  phase  steels  are  in  agreement  with 
the  theory  of  Mileiko  (5)  and  Garmong  and  Thompson  (6).  Since  this 
theory  (5,6)  has  been  developed  for  mechanical  properties  of  fiber 
composites  of  two  ductile  phases  tested  in  tension  parallel  to  the  fiber 
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axis,  it  will  have  only  a  limited  applicability  because  it  (5,6)  does  not 
consider  the  morphology  of  the  phases,  which  can  effect  the  stress- 
strain  behavior,  as  will  be  shown  later. 


Tomota  et  al.  (7)  have  calculated  the  flow-stresses  of  alloys  con¬ 
sisting  of  two  ductile  phases,  a-y  Fe-Cr-Ni  alloys  on  the  basis  of  a 
continuum  model  where  the  internal  stress  produced  by  the  inhomo¬ 
geneous  distribution  of  plastic  strain  is  taken  into  account.  They  (7) 
have  found  that  the  calculated  curve  of  the  two-phase  alloy  was  al¬ 
most  identical  with  that  constructed  by  the  simple  equal  strain  model 
but  the  authors  caution  that  there  are  circumstances  where  the  mater¬ 
ial  would  not  behave  as  if  constant  strain  existed. 

Fischmeister  et  al.  (8),  Karlsson  and  Sundstrom  (9),  and 
Fischmeister  and  Karlsson  (10)  have  recently  applied  the  FEM  (11)  to 
calculate  the  stress-strain  curve  of  the  soft  matrix-hard  particle, 
ferrite-martensite  structures  with  30  volume  percent  of  martensite. 
In  a  previous  investigation  (12),  FEM  has  been  applied  to  calculate 
the  stress-strain  curve  of  a  soft  particle-hard  matrix,  a-0  Ti-Mn  alloy 
with  16.6  vol.%  a,  and  it  was  shown  that  the  calculated  curve  lay 
below  the  experimental  stress -strain  curve.  Fig.  7.  This  difference 
was  attributed  to  two  factors.  The  first  was  the  fact  that  the  stress- 
strain  curves  of  "single”  phase  a  and  p  phases  were  used  for  the 
calculations  (12).  Since  the  grain  sizes  of  both  the  "single  phase"  a 
and  8  was  larger  than  the  grain  sizes  of  these  phases  in  the  duplex 
alloy  (13),  the  calculated  curve  might  have  been  lower  than  the 
experimental  curve.  The  second  factor  was  that  the  presence  of  the 
interface  phase  (14-16)  might  have  contributed  to  the  strength  of  the 
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duplex  alloy  (17)  thereby  raising  the  experimental  curve  above  the 
calculated  curve,  which  did  not  consider  this  structure.  In  the 
present  calculations  the  first  factor  has  been  considered  and  the 
stress-strain  curve  has  been  recalculated.  The  second  factor  is  not 
considered,  since  it  is  expected  that  the  interface  phase  would  make 
only  a  small  contribution  to  the  mechanical  properties  of  a  -p  Ti-Mn 
alloys.  Other  factors,  namely  the  effect  of  particle  size  and  matrix 
which  were  not  considered  in  the  previous  analysis  (12),  have  been 
taken  into  account  in  the  present  analysis.  In  addition  the  calcu¬ 
lations  were  made  not  only  for  the  16.3  vol.%  a  alloy  but  also  for  the 
80,  63  and  41  vol.  %  a,  and  all  the  calculations  have  been  compared 
with  their  corresponding  experimental  curves. 

PROCEDURE 

The  procedure  used  in  the  present  calculations  is  essentially 
similar  to  the  one  used  in  the  previous  calculations  (12)  and  hence 
will  be  described  only  briefly  here.  However,  there  are  some 
modifications  and  a  check  which  will  be  described.  NASTRAN  (18) 
computer  program,  in  particular  rigid  format  No.  6  called  "piecewise 
linear  analysis",  has  been  used  for  the  solution  of  the  elasto-plastic 
problem.  An  IBM  360/65  computer  has  been  used.  For  all  the 
calculations  a  uniform  mesh  of  392  triangular  two  dimensional  (plane 
stress)  plate  elements.  Fig.  1  was  used,  and  the  volume  fraction, 
particle  size  and  shape  could  be  varied  by  designating  e<.ch  triangle 
either  a  or  p,  Figs.  1  to  6.  The  node  J  in  Figs.  1  to  6  is 
considered  fixed  and  all  the  other  nodes  along  the  line  ^  can  move 
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only  in  the  x  direction.  All  the  other  nodes  can  move  either  in  the  x 
or  y  directions. 

In  the  previous  calculation  (12)  the  nodes  along  the  line  CD 

moved  independently  m  tne  x  and  y  directions.  Excessive  differences 
in  the  y  displacements  in  this  edge  CD  were  found  (12)  beyond  the 
stress  level  758  MPa  (110  KSi),  which  is  not  expected  in  the  real 
specimens.  Thus  all  the  calculations  had  to  be  restricted  to  758  MPa 
(12).  In  the  present  analysis  this  problem  was  overcome  through  a 
multipoint  constraint  such  that  the  nodes  along  the  edge  CD  have  the 
same  y  displacement  at  any  stage  of  deformation  whereas  all  the  nodes 
other  than  the  nodes  in  the  lines  AB  and  CD  still  can  move  freely 
either  in  x  or  y  directions.  This  multipoint  constraint  was  not  incor¬ 
porated  in  the  calculations  of  the  matrix  effect.  Fig.  4,  where  the 
previous  procedure  (12)  was  followed. 

The  input  stress-strain  curves  for  a  and  3  phases  that  were 
used  for  all  the  calculations,  except  for  the  calculations  of  the  matrix 
effect.  Fig.  4  and  9,  are  shown  in  Fig.  7.  The  a  stress-strain  curve 
which  corresponds  to  an  a  grain  size,  D,  of  3ym*  has  been  derived 
from  the  stress-strain  curves  of  larger  grain  sizes,  by  extrapolating 
the  flow  stress  vs.  D"^2  plots  for  specific  strains.  It  is  to  be  noted 
that  the  flow  stress  dependency  on  grain  size  of  or  was  found  to  be 
small  and  is  close  to  that  found  by  Jones  and  Conrad  (19)  for  Iodide 
of  titanium,  even  though  the  authors  a  titanium  had  an  oxygen 

*The  a  grain  size  of  the  alloys,  the  stress-strain  curves  of  which  are 
shown  in  Fig.  12,  varied  from  3  to  6  pm.  This  variation  in  the  a 

grain  size  does  not  significantly  effect  the  calculated  curves  which 

are  also  shown  in  Fig.  12. 


content  of  800  p.p.m.  The  flow  stress  dependency  on  grain  size  of  p 
titanium  was  found  to  be  negligibly  small.  However,  in  the  present 
calculations  the  stress-strain  curve  of  p,  Fig.  7,  without  any  a  was 
used,  unlike  the  previous  calculations  (12)  where  the  p  had  about  4 
vol.  %  a.  In  brief,  the  a  and  p  curves  used  in  the  present  calcula¬ 
tions,  Fig.  7,  are  slightly  higher  than  the  curves  used  in  the  pre¬ 
vious  calculations  (12).  For  example,  at  a  total  strain  of  0.015,  the 
flow  stress  difference  for  a  curves  is  59  MPa  (8.5  KSi)  and  the  flow 
stress  difference  for  p  curves  is  31  MPa  (4.5  KSi);  while  the  differ¬ 
ences  in  the  elastic  part  of  the  curves  are  negligibly  small.  As 
mentioned  before,  the  a  and  p  curves  of  the  previous  calculations 
(12)  were  used  in  the  calculations  of  matrix  effect,  Figs.  4  and  9. 

The  load  was  applied  uniformly  at  the  nodes  of  the  edge  CD  and 
was  increased  by  predetermined  increments  where  the  first  loading 
condition  is  such  that  all  the  stresses  would  be  in  the  elastic  region. 
For  plotting  the  stress-strain  curves.  Figs.  7,  8,  10  and  12,  the 
stresses  are  known  for  each  loading  condition,  and  the  strains  are 
calculated  from  the  common  displacement  of  the  nodes  of  the  line  CD 
and  the  initial  length  AC,  Figs.  I,  2,  3,  5  and  6.  For  calculating 
the  stress-strain  curve  of  the  matrix  effect  (Fig.  9),  the  stresses  are 
known  for  each  loading  condition  and  the  strain  was  calculated  from 
the  weighted  average  of  the  displacements  of  the  nodes  of  the  line  EF 
and  the  initial  length  AE,  Fig.  4. 

To  check  the  reproducibility  of  the  method  of  calculations,  with 
the  present  modifications  and  procedure,  and  also  to  find  out  the  limit 
of  the  strain  up  to  which  calculations  could  be  performed  with  mini¬ 
mum  error,  the  whole  mesh  in  Fig.  1  was  made  a  and  the  calculations 


were  performed.  It  was  found  that  the  calculated  a  curve  closely 
followed  the  input  a  curve  upto  a  strain  of  0.05,  after  which  the  cal¬ 
culated  a  curve  rose  sharply  above  the  input  a  curve,  not  shown 
here.  Hence,  all  the  calculations  were  restricted  to  a  total  strain 
0.05  or  less  in  the  a-p  alloys  considered.  It  should  be  noted  that 
even  though  the  average  strains  in  the  a-p  alloys.  Figs.  7  to  10  and 
12,  are  lower  than  0.05  strain,  there  are  instances,  in  particular  for 
stress  levels  896  MPa  (130  KSi)  and  965  MPa  (140  KSi)  of  the  config¬ 
uration  of  Fig.  1  and,  for  the  stress  level  965  MPa  (140  KSi)  of  the 
configuration  of  Fig.  2,  for  which  the  average  strain  in  the  a  phase 
is  slightly  higher  than  0.05.  Since  the  volume  percent  of  a  in  these 
configurations,  Figs.  1  and  2,  is  only  16.3%,  the  errors  arising  from 
this  strain  would  have  only  a  small  effect  on  the  calculated  a-p 
stress-strain  curves.  Fig.  8. 

To  ascertain  the  effect  of  particle  size  on  the  calculated  or-e 
curve,  the  particle  size  was  varied  as  shown  in  Figs.  1  to  3,  while 
keeping  the  16.3  vol.  %  a  constant. 

The  mesh  used  to  calculate  the  effect  of  matrix  is  shown  in  Fig. 
4.  For  this  calculation  only  the  area  of  the  mesh  ABFE  was  consid¬ 
ered  and  the  area  EFDC  was  made  as  the  stronger  phase,  p,  so  that 
it  served  as  approximately  a  rigid  body,  with  minimum  differences  in 
the  displacements  of  the  nodes  along  the  edge  CD.  This  had  to  be 
done  because  the  multipoint  constraint  described  above  was  not  incor¬ 
porated  for  this  calculation,  as  mentioned  above.  Within  the  area 
ABFE,  the  volume  fractions  of  a  and  p  are  0.50.  In  one  calculation 
the  shaded  phase.  Fig.  4,  was  made  a,  the  particles  and  the  un¬ 
shaded  phase  p  as  the  matrix.  In  the  comparison  calculation  the 
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designations  were  reversed.  In  this  process  of  interchanging  the  two 
phases,  nothing  else  was  changed,  and  hence  the  difference  between 
the  two  calculated  curves,  Fig.  9,  gives  the  effect  of  the  matrix. 

Calculations  were  also  made  for  the  volume  percents  of  a  of  41, 
63  and  80  and  for  the  configurations  of  the  particle  si2es  similar  to 
Figs.  1  and  3.  However,  the  meshes  used  for  the  calculations  for  the 
41  and  63  volume  percents  a  are  not  shown  here.  When  the  volume 
percent  of  or  was  80,  the  p  was  in  the  form  of  particles  and  the 
configurations  of  p  for  this  a  -  p  alloy  are  shown  in  Figs.  5  and  6. 

The  calculations,  made  with  the  minor  phase  in  the  form  of 
triangles,  as  in  Fig.  3,  are  compared  with  the  experimentally  deter¬ 
mined  curves  for  all  four  a  -  p  alloys  in  Fig.  12.  All  the  four  a-p 
alloys  were  annealed  for  200  hrs.  at  700°C  (973°K),  and  tensile  tests 
were  conducted  on  a  Tinius  Olsen  testing  machine  with  2.54  cm  (1 
inch)  standard  gauge  length  specimens. 


RESULTS  AND  DISCUSSION 
A.  Comparison  with  the  Previous  Work  (12): 


The  calculated  curve,  obtained  with  the  mesh  of  Fig.  1  for  the 
16.3  vol.  %  a  alloy  is  shown  in  Fig.  7.  This  stress-strain  curve  lay 
above  the  calculated  curve  of  previous  work  (12),  Fig.  7,  for  the 
same  volume  percent  of  a.  The  flow  stress  difference  between  the 
two  curves  at  a  true  strain  of  0.015  is  79  MPa  (11.5  KSl).  Part  of 
this  difference  obviously  is  due  to  the  difference  in  input  stress- 


strain  curves  of  the  a  and  p  phases  for  the  two  calculations.  Part 
of  this  difference  could  also  be  due  to  the  size  and  shape  of  the  or 
particles  used  in  the  mesh  of  Fig.  1  and  in  the  mesh  of  the  previous 


work  (12),  even  though  the  interface  area/unit  volume  is  not  signi¬ 
ficantly  different. 

Even  though  the  calculated  curve  of  the  present  work  is  above 
the  calculated  curve  of  the  previous  work  (12),  it  is,  nevertheless 
considerably  below  the  experimental  curve,  shown  in  Fig.  7.  The 
flow  stress  difference  between  these  two  curves  at  a  strain  of  0.015 
is  about  104  MPa  (15KS1).  If  one  recognizes  that  in  the  mesh  of 
Fig.  1.  AB  represents  the  diameter  of  the  specimen,  6.35  mm,  then 
the  particle  size  in  the  mesh  is  very  much  larger  than  the  particle 
size  of  the  test  specimen.  Thus  it  appeared  desirable  to  ascertain 
whether  the  discrepancy  between  calculated  and  experimental  results 
was  related  to  particle  size. 

B.  Effect  of  Particle  Size: 

The  meshes  used  to  calculate  the  effect  of  particle  size  for  a 
constant  16.3  vol.  %  a  are  shown  in  Figs.  1  to  3  and  the  corres¬ 
ponding  calculated  stress-strain  curves  are  shown  in  Fig.  8.  It  can 
be  seen  that,  as  the  particle  size  gets  smaller,  the  calculated  stress- 
strain  curve  tends  to  be  higher,  Fig.  8.  The  variation  in  interface 
area/unit  volume  is  3.5  times  when  the  particle  size  varied  from 
coarse,  Fig.  1  to  fine,  Fig.  3.  The  flow  stress,  for  a  given  strain, 
for  the  fine  a  particle  configuration  is  higher  than  that  for  the  coarse 
a  particles,  because  the  strain  differences  between  the  a  and  p 
phases  are  smaller  for  the  fine  than  for  the  coarse  particles.  The 
smaller  strain  difference  means  that  the  p  has  also  undergone  consid¬ 
erable  larger  strains  which,  in  turn,  means  a  higher  flow  stress, 
because  p  is  stronger  than  or. 


The  difference  in  the  flow  stress  between  the  calculated  curves. 
Fig.  8,  starts  to  be  different,  almost  from  the  point  of  proportional 
limit,  414  MPa  (60KSi).  Thereupon  the  difference  between  the  flow 
stresses,  for  specific  strains,  tends  to  be  increasingly  different  up  to 
a  strain  of  0.0125,  at  which  point  the  maximum  difference  occurs. 
This  maximum  difference  between  the  coarse  and  fine  particles,  upper 
and  lower  curves  in  Fig.  8,  is  about  69  MPa  (10  KSi)  and  this  differ¬ 
ence  remains  fairly  constant  up  to  a  strain  of  0.0262.  Comparison  of 
the  calculated  curves.  Fig.  8  for  coarse.  Fig.  1,  and  medium  size. 
Fig.  2,  a  particles  indicates  that  the  difference  in  the  flow  stress 
tends  to  be  reduced  beyond  the  strain  0.0262  and  the  flow  stresses 
are  very  close  at  a  strain  0.035. 

An  attempt  has  been  made  to  determine  what  size  of  the  a  parti¬ 
cles  in  the  mesh  needs  to  be  used  so  that  the  calculated  curve  closely 
approaches  the  experimental  curve.  This  was  done  by  extrapolating 
the  flow  stress  vs.  the  particle  size**  for  specific  total  strains  to  the 
flow  stress  of  the  experimental  curve,  Fig.  7.  It  was  found  that  this 
size  depends  on  the  magnitude  of  the  total  strain.  The  extrapolated 
particle  size  at  a  total  strain  of  0.0075  was  30  pm  and  150  pm  at  a 
total  strain  of  0.025,  these  sizes  are,  respectively,  1/9  and  3/5  of  the 
particle  size  of  Fig.  3.  At  intermediate  total  strains,  the  particle 
sizes  were  intermediate  also.  It  is  to  be  noted  that  the  size  of  the  a 
particles  in  the  actual  material  was  4.4  pm.  Therefore,  these  results 
suggest  that  at  strains  above  0.025  the  particle  size  (or  the  interface 

**  If  the  line  AB  in  Figs.  1  to  3  is  considered  as  the  diameter,  of 
the  specimen,  then  the  equivalent  particle  sizes  in  these  figures 
are  906,  453  and  259  pm,  respectively. 
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area/unit  volume)  does  not  have  much  influence  on  the  flow  stress, 
whereas  at  lower  strains  particle  size  does  influence  the  flow  stress. 
Some  support  for  this  line  of  reasoning  comes  from  the  work  of  Mar¬ 
golin  and  Stanescu  (20).  These  workers  showed  for  polycrystailine 
material,  that  the  difference  in  the  strength  levels  between  the  grain 
boundary  zone  and  grain  interior  is  largest  at  low  plastic  strains  and 
decreased  with  increasing  plastic  strains  until  at  2%  the  two  regions 
had  the  same  flow  stress. 

As  will  be  shown  later.  Fig.  12,  the  calculated  curve  with  fine  a 
particles  is  close  to  the  experimental  curve.  Therefore,  these  calcula¬ 
tions  underline  the  need  to  use  fine  a  particles,  or  a  large  area  of 
the  actual  microstructure  to  represent  the  a  particles  in  the  actual 
a  -  0  alloys,  a  need  not  expressed  in  previous  investigations 
(8-10,12). 

C.  Effect  of  Matrix: 

The  calculated  stress-strain  curves  are  shown  in  Fig.  9.  It  can 
be  seen  that  the  curve  with  0  as  matrix  lies  above  the  the  curve  with 
a  as  matrix.  This  behavior  is  expected,  since  0  is  stronger  than  a. 
Since  coarse  particles  were  used  for  these  calculations  it  is  to  be 
expected  that  both  curves  would  be  raised  with  finer  particle  sizes. 
However,  as  can  be  seen  from  a  comparison  of  Figs.  8  and  10,  the  0 
matrix  curve  would  be  expected  to  be  raised  more  than  the  a  matrix 
curve.  Comparison  of  the  two  curves  in  Fig.  9  indicates  that  the 
flow  stress  differences,  at  a  given  true  strain,  starts  to  be  signi¬ 
ficantly  different  only  after  about  0.2%  plastic  strain,  unlike  the 
effect  of  particle  size  where  the  flow  stress  differences  started  to  be 
significant  from  the  proportional  limit  onward.  Fig.  8. 
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D.  Effect  of  Volume  Fraction: 


Calculations  for  coarse  and  fine  particle  sizes,  comparable  to 
Figs.  1  and  3  were  also  made  for  41,  63,  and  80  vol.  %  a  alloys.  The 
calculated  stress -strain  curves  for  the  80  vol.  %  a  alloy  are 
shown  in  Fig.  10.  Figure  10  reveals  that  the  difference  between  the 
calculated  curves  starts  to  be  significant  right  from  the  proportional 
limit,  as  was  found  earlier  for  the  16.3  vol.%  a  alloy,  Fig.  8,  and 
reaches  a  maximum  which  is  48  MPa  (7  KS1)  at  a  strain  0.0232. 

The  behavior  of  the  calculated  curves  for  the  coarse  and  fine 
second  phase  particles  for  the  63  and  41  vol.%  a  alloys  were  in  be¬ 
tween  that  of  the  80,  Fig.  10,  and  16.3,  Fig.  8.  vol.%  a  alloys  and 
hence  are  not  shown  here.  However,  the  0.2%  and  0.4%  offset  YS 
differences  between  the  coarse  and  fine  second  phase  particle  cal¬ 
culated  curves,  for  all  the  four  a-p  alloys  are  given  in  Fig.  11.  It 
indicates  that  the  particle  size  effects  are  more  important  for  the  YS 
determination  in  the  range  of  35  to  90  volume  percent  p. 

E.  Comparison  of  Calculated  and  Experimental  Curves: 

The  calculated  and  experimental  stress-strain  curves  for  the  four 
a  -  p  alloys  are  presented  in  Fig .  12 .  The  calculated  curves  are  for 
fine  second  phase  particles.  It  can  be  seen  that  for  all  the  alloys, 
the  FEM  calculated  curves  predict  deviation  from  elastic  behavior  at  a 
lower  stress  than  the  experimental  curves.  In  addition,  the  calcula¬ 
ted  flow  stresses  at  strains,  just  after  the  proportional  limits  of  the 
experimental  curves,  are  lower  than  the  experimental  flow  stresses. 
For  example,  for  the  alloys  2,  3,  4  and  5  at  the  strains  0.0039, 
0.005,  0.0063  and  0.0101,  the  differences  are  48  MPa  (7  KSi),  79  MPa 
(11.5  KSi),  69  MPa  (10  KSi)  and  48  MPa  (7  KSi),  respectively. 

-44- 


These  are  the  maximum  differences  between  the  flow  stresses  of  the 
calculated  and  experimental  curves  and  at  higher  strain  levels,  the 
differences  are  reduced.  These  calculations  are  consistent  with  the 
FEM  calculation  of  Fischmeister  and  Karlsson  (10)  where  it  was  ob¬ 
served  that,  for  a  ferrite  +  martensite  steel,  the  maximum  difference 
between  the  calculated  and  experimental  results  occur  at  the  very 
beginning  of  the  plastic  deformation. 

As  mentioned  earlier,  the  particle  sizes  of  Figs.  3  and  6  are 
very  much  coarser  than  the  actual  particle  sizes  of  the  alloys.  Conse¬ 
quently,  if  still  finer  particle  sizes  were  used,  it  would  be  reasonable 
to  expect  that  the  calculated  curves  would  approach  the  experimental 
curves  more  closely  still.  Further  support  for  this  view  comes  from 
the  fact  that  of  all  the  four  alloys,  mentioned  above,  the  maximum 
difference  of  79  MPa  (11.5  KSi)  between  calculated  and  experimental 
curves  occurred  for  alloy  3,  which  also  was  found  experimentally  to 
have  the  highest  interface  area/unit  volume.  The  differences  in  the 
calculated  and  experimental  curves,  in  general,  are  very  small  at 
higher  strains.  Fig.  12. 

Fig.  12  shows  that  at  the  larger  plastic  strains,  the  calculated 
curve  lies  above  the  experimental  curve  for  alloy  2 ,  whereas  for 
alloys  3,  4  and  5  alloys  the  calculated  curves  lie  below  the  experimen¬ 
tal  curves.  As  earlier  results  on  the  matrix  effect  for  a  50  vol.%  a 
revealed,  when  0  was  the  matrix,  the  calculated  stress-strain  curve 
was  higher  than  when  a  was  the  matrix.  The  p  particles  arranged  in 
the  mesh  of  Fig.  6  for  alloy  2  tend  to  be  more  continuous  than  the  p 
particles  in  the  actual  material.  Thus,  the  calculated  curve  would  be 
expected  to  lie  above  the  experimental  curve  at  the  larger  strains,  as 


shown  in  Fig.  12.  For  the  alloys,  3,  4  and  5,  the  a  particle  arrange¬ 
ment  tend  to  be  more  continuous  than  a  in  the  actual  materials  and 
hence  lower  calculated  curves  than  the  experimental  curves  would  be 
expected.  This  is  also  seen  in  Fig.  12.  With  the  exceptions  noted, 
it  appears  that  the  FEM  satisfactorily  predicts  the  stress-strain  behav¬ 
ior  of  the  two  ductile  phase  alloys,  from  the  stress-strain  behavior  of 
the  component  phases  and  volume  fractions. 

It  should  be  noted,  however,  that  in  addition  to  the  0.05  strain 
limitation  mentioned  earlier,  the  present  use  of  the  FEM  has  not  con¬ 
sidered  the  anisotropy  of  behavior  of  the  alpha  and  beta  phases, 
either  elastically  or  plastically.  No  distinction  has  been  made  between 
slip  behavior  in  equiaxed  and  Widmanstatten  alpha.  Furthermore  the 
presumption  of  compatibility  at  a-8  interfaces  rules  out  the  possibility 
of  a  breakdown  of  the  compatibility  which  would  lead  either  to  void 
formation  or  cracking  which  can  occur  at  very  low  strains  (21). 

CONCLUSIONS 

1  By  using  NASTRAN  (18)  computer  program  and  stress-strain 
curves  of  the  a  and  £  phases  of  Ti-Mn  alloys,  the  Finite  Element 
Method  (FEM)  has  been  successfully  employed  to  calculate  the 
effect  of  particle  size,  matrix  and  volume  fraction  on  the  stress- 
strain  relations  of  two  phase  of-0  titanium  alloys. 

2)  Calculations  of  the  stress-strain  curve  of  a  16  volume  percent  a 
alloy  were  carried  out  for  three  particle  sizes.  As  the  particle 
size  decreased,  the  calculated  curve  approached  the  experimental 
curve  and  was  close  to  it  for  the  smallest  particle  size  used,  260 
pm.  Calculations  showed  that  the  particle  size  which  would 
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produce  the  experimental  curve  depended  upon  the  strain,  and 
the  smallest  particle  sizes  were  required  for  the  smallest  strain. 
In  general  deviation  between  calculated  and  experimental  curves 
for  the  80,  63,  41  and  16.3  volume  percent  a  alloys  were  largest 
at  the  onset  of  plastic  deformation  and  diminished  with  increasing 
strain.  Additional  calculations  showed  that  the  particle  size 
effect  was  strongest  in  the  region  35-90  volume  percent  p. 

3)  The  matrix  effect  was  determined  by  using  a  mesh  in  which  the 
a  and  |3  phases  were  present  in  equal  amounts.  The  properties 
of  the  matrix  phase  were  alternately  assigned  to  a  and  6.  It 
was  found,  as  anticipated,  that  when  the  stronger  phase,  p,  was 
matrix,  the  stress-strain  curve  reached  higher  stresses.  Be¬ 
cause  the  particle  size  used  was  coarse,  it  is  quite  likely  that 
the  strengthening  effect  of  the  p  matrix  phase  was  not  fully  rea¬ 
lized.  The  small  differences  between  calculated  and  experimental 
curves  were  found  to  be  consistent  with  the  "degree  of  ma- 
tricity"  of  the  particle  phase  as  well  as  the  particle  size. 
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FIGURES 

Fig.  1.  Mesh  1,  Coarse  a  particles  (16.3  vol.  %),  shaded,  in  a  p 
matrix,  unshaded. 

Fig.  2.  Mesh  2,  Medium  size  a  particles  (16.3  vol.  %),  shaded,  in 
the  p  matrix,  unshaded. 

Fig.  3.  Mesh  3,  Fine  a  particles  (16.3  vol.  %),  shaded,  in  the  p 
matrix,  unshaded. 

Fig.  4.  Mesh  4,  Particles  Phase  shaded  and  matrix  phase  unshaded. 

Particles  and  Matrix  in  equal  amounts  within  the  area  of  the 
mesh  ABFE. 

Fig.  5.  Mesh  5,  Coarse  p  particles  (20  vol.  %),  shaded,  in  the  a 
matrix,  unshaded. 

Fig.  6.  Mesh  6,  Fine  p  particles  (20  vol.  %),  shaded,  in  the  a 
matrix  unshaded. 

Fig.  7.  The  stress-strain  curves  of  the  a  and  p  phases  used 

for  the  FEM  calculations,  calculated  curve  for  a  16.3  vol.  % 
a  alloy  with  coarse  a  particles  with  Mesh  1  of  Fig.  1,  FEM 
calculated  curve  for  a  16.6  vol.  %  a  alloy  from  the  previous 
work  (12)  and,  experimentally  determined  stress-strain 
curve  for  a  similar  (~  17)  vol.  %  a  alloy. 

Fig.  8.  FEM  calculated  stress-strain  curves  for  a  16.3  vol.  %  a 

alloy  with  three  different  a  particle  sizes.  Figs.  1  to  3.  . 

Fig.  9.  FEM  calculated  stress-strain  curves  for  a  50  vol.  %  <*-p 
alloy,  Fig.  4,  one  with  a  and  p  alternately  as  matrix. 

Fig.  10.  FEM  calculated  stress- strain  curves  for  an  80  vol.  %  a 

alloy  with  two  different  particle  sizes  of  p,  Figs.  5  and  6. 

Fig.  11.  0.2  and  0.4%  offset  yield  strength  differences  between  the 

coarse  and  fine  second  phase  particles  vs.  the  volume  per¬ 
cent  of  p. 

Fig.  12.  Comparison  of  the  FEM  calculated  stress-strain  curves,  with 
fine  second  phase  particles,  with  corresponding  experimental 
stress-strain  curves. 
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FIGURE  1 


FIGURE  2 
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